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Q | In entangled polymer systems, there are several characteristic time scales, such as the entangle- 

, ment time and the disengagement time. In molecular simulations, the longest relaxation time (the 

ry^ ■ disengagement time) can be determined by the mean square displacement (MSD) of a segment or 

by the shear relaxation modulus. Here, we propose the relative fluctuation analysis method, which 
is originally developed for characterizing large fluctuations, to determine the longest relaxation time 
from the center of mass trajectories of polymer chains (the time-averaged MSDs). Applying the 
method to simulation data of entangled polymers (by the slip-spring model and the simple rep- 
tation model), we provide a clear evidence that the longest relaxation time is estimated as the 
C/j , crossover time in the relative fluctuations. 



INTRODUCTION 



Polymer melts and solutions exhibit interesting dynamical behavior. In particular, if the degree of polymerization is 
large, characteristic dynamical behavior due to the entanglement effect can be observedp]]. For example, viscoelasticity 
, with a very long relaxation time (which strongly depends on the molecular weight) and anomalous diffusion of a 
segment occur as a result of the entanglement effect. Dynamics of entangled polymers have been studied extensively 
by various methods (lH3|, including theoretical modeling and simulations. 

There are several characteristic time scales in entangled polymer systems, such as the entanglement time r e , the 
Rouse time tr, and the disengagement time (the longest relaxation time) Td- In each characteristic time scale, mean 
square displacements (MSDs) show different scalings. For example, the Doi-Edwards tube model gives the following 
time dependence of the MSD of a single segment [l[: 



[t ,/s (t<r.) 

[t (Td < t), 



where r.;(i) is the position of the i-th segment at time t and (...) represents the ensemble average. Anomalous 
diffusion of a single segment in the Doi-Edwards tube model originates from the following two mechanisms. One is 
a dynamical property of a single Rouse chain. The MSD of a segment of a Rouse chain shows subdiffusive behavior, 
([ r i(t) — fi(0)] 2 ) oc t 1 / 2 , at the short time scale (t < i7j)[l|. The other is a constraint by a tube (which is formed by 
surrounding chains). Namely, each segment can move only along the tube, which has a fractal geometry similar to a 
random walk trajectory, thereby causing a subdiffusive transport characterized by the exponent 1/2 for r e < t < Td- 
The subdiffusive exponent 1/4 for r e < t < tr results from the combination of these two mechanisms. Thus, using 
the MSD of a single segment, one can estimate r e , tr and 7^0-01- 

It should be noticed that numerical prefactors in Eq. ([1]) (which are not explicitly shown) depend on the index 
i if the target segment is close to chain ends[6(. (And if the strong z-dependence is observed, validity of Eq. ([]} is 
no longer guaranteed^ Although this does not cause serious problems for relatively fine-scale models (such as the 
Kremer-Grest model[3-@), it may be crucial for highly coarse-grained models (such as the RaPiD (responsive particle 
dynamics) model[7, 8]). In some coarse-grained models, resolutions of the models are not sufficient to resolve the 
MSDs of segments. Therefore, for such coarse-grained models, it is physically reasonable to employ the MSD of the 
center of mass (CM) of a chain, instead of that of a single segment. In contrast to the MSD of a single segment, the 
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MSD of the CM is given by 

U (t< r e ) 

([rcM(i)-^M(0)]V It 1 ' 2 (r e <t<r R ) (2) 
I* (TR<t), 

where rcM(i) is the position of the CM. We note that the scaling exponent in the MSD of the CM does not change 
around t fts t<j. Roughly speaking, this is because the constraint by a tube apparently disappears by taking the 
average with respect to all segments. It seems to be difficult to extract an information of from the MSD of the CM 
straightforwardly. 

For highly coarse-grained models, thus, usually other physical quantities such as the shear relaxation modulus 
and the end-to-end vector relaxation function have been employed to determine the long time relaxation behavior. 
However, it is not always easy to determine the long time relaxation behavior accurately. For example, to obtain 
the accurate shear relaxation modulus data, we need to compute the autocorrelation functions of stress tensors for 
large number of statistically independent samples. Such calculation becomes demanding for interacting many chain 
models. The end-to-end vector relaxation function can be computed more accurately than the relaxation modulus, 
but some highly coarse-grained models such as the RaPiD model do not provide the information of end-to-end vectors. 
Therefore, alternative methods to determine the long time behavior will be useful for some simulation models, and 
thus are worth developing. A simple but fascinating idea is to utilize some tools or methods developed in the fields 
other than polymer physics to analyze simulation data of entangled polymers. 

In this work, we employ a method which is originally developed in Levy statistics 9] , to analyze the long time 
relaxation behavior of entangled pol ym ers. Recently, various theoretical tools have been developed for the continuous 
time random walk (CTRW) model 10], which is a typical trap model. An interesting recent finding is that the 
distribution of the time-average MSD (TAMSD) can be utilized to characterize some dynamical properties [Ill - fill . 
The TAMSD for a position r is defined as 

1 pt~A 



^ (A; i) = __ ^ ^>(t' + A)-r(i')] 2 , (3) 
with A and t being the time difference and the observation time. The relative fluctuation (RF) of the TAMSD 



R(t; A) = (4) 
(S 2 (A;t)) 

is useful for characterizing non-Gaussian fluctuations and ergodicity breaking in diffusion processes [ill [l3l-[T5|. For 
the CTRW trajectories, by analyzing the i-dependence of the RF, we can determine some important information of 
the trapping-time distribution, such as the power-law exponent and relaxation timefl3l. Il5j|. 

Although the CTRW is not equivalent to models for entangled polymers, the RFs in polymer systems may also 
characterize some long time relaxation behavior. If the RF analysis (RFA) works well for entangled polymers, it will 
serve as a new and complementary tool to study the long time relaxation behavior. As far as the authors know, such 
an analysis has never been performed for entangled polymer systems. In this work we perform the RFA for the CM 
trajectory data of entangled polymer models and show that the information of the long time relaxation behavior (such 
as the longest relaxation time Td) can be observed. Therefore the RFA can be utilized as a new and complementary 
method. 



II. SIMULATION RESULTS 



In this section, we perform simulations for entangled polymers and analyze MSDs of CMs. So far, various simulation 
models have been proposed and utilized to study entangled polymers d, 0, B Il7l - [2"2l |. In this work, we employ two 
simple models among them. One is the slip-spring model[21J, |22| and the other is the simple Doi-Edwards type 
reptation model[23[. We calculate the RFs of the TAMSDs, Eq. flU), for the CM trajectories (t"cm(*)), and show that 
the RFs show clear crossovers around t ~ Td- 



A. Slip-Spring Model 

First we perform simulations of the slip-spring model[2l|, [22[. The slip-spring model is a variant of slip- link 
models [17H20(, which reproduce diffusive and rheological properties well. In the slip-spring model, a polymer chain is 
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described as an ideal chain and the entanglements are mimicked by so-called slip-springs. One end of a slip-spring is 
anchored in space while another end can slip along the chain. Moreover, slip-springs are dynamically reconstructed 
at chain ends. This model can reasonably reproduce linear and nonlinear rheological properties, while it lacks some 
mechanisms such as the constraint release (CR) or the convective constraint release (CCR)@, Q. We perform simula- 
tions for different numbers of beads (polymerization indices) TV. Other parameters are set to be the standard values; 
the average number of segments between slip-springs Nq — 4, the slip-spring strength N s = 0.5, the friction coefficient 
of a slip-spring £ s /£ = 0.1 (with £ being the friction coefficient of a segment). As elementary length and time scales 
of this model, we use the segment size b and the characteristic time of a segment tq = £6 2 /ksT. (ks and T are 
the Boltzmann constant and the temperature, respectively.) For each TV value, simulations are performed for 1024 
polymer chains sampled from an equilibrium distribution. (Details of simulation model and algorithm used in this 
work are described in Ref . [22| ) 

Figure ED shows the ensemble-averaged MSDs (EAMSDs) of CMs for TV = 10,20,40,80, and 160. (The CM is 
determined as the average position of the beads which compose a polymer chain.) For large TV cases, the t dependence 
of the MSDs is well described by Eq. (jS). Thus information of two characteristic times, r e and tr, can be obtained. 
However, as we mentioned in Section HJ we cannot obtain any information about Td from the EAMSDs for CMs. Here 
we perform the RFA of the TAMSDs for CMs. As will be shown in Appendix [A] it is difficult to obtain the RF 
accurately if A is too small. Also, if A is too large, the RF shows only Gaussian fluctuation. Therefore, the time 
difference A should be chosen to be some moderate values. Here, we employ A/tq = 10 for the current slip-spring 
simulations. 

Figure |2] shows the RF for different TV. For large TV cases (TV > 40), the RF is fitted well by the following form 
except for a very short time region: 

where a is a power-law exponent (0 < a < 0.5 ) and r c is the crossover time. The crossover time, r c , is much 
longer than the crossover time for the EAMSD and strongly depends on TV. We cannot determine the crossover time 
accurately for small TV (TV < 20, Fig.[2£a)). For TV > 40 (Fig.^b)), we can determine the crossover time from the RF 
data easily. The crossover time r c strongly depends on TV. (The power-law exponent a also depends on TV. Fitting 
gives a = 0.31, 0.24 and 0.19 for TV — 40, 80 and 160, respectively. We will discuss this in Section Hill ) Such behavior 
of t c seems to be consistent with the behavior of the longest relaxation time Td determined from the shear relaxation 
modulus [22| . Actually, the TV dependence of t c obtained from the RFA is quite similar to that of Td (see Fig. For 
TV > 40, Td and t c show the similar power-law type TV dependencies (r<j oc TV 3 - 48 and r c cx TV 3 ' 51 ). 

For TV < 20, the TV dependence of Td is close to one of unentangled polymers (t^ oc TV 2 ) and in this region r c is not 
determined accurately (see Fig. Ufa)). Possibly, this implies that we cannot extract the information about the long 
time relaxation behavior from the RFA if Td is comparable to or less than A. Similar trend is found in Appendix [A] 
Although behavior of the RF of the TAMSDs depends on parameters such as TV and A, we have qualitatively the 
same result for other cases as long as TV is sufficiently large. Thus, in the slip-spring model, we conclude that r c is 
essentially the same as Td apart from the numerical constant. 

B. Discrete Reptation Model 

To confirm whether the relation between the cutoff time and the longest relaxation time is a general feature in 
entangled polymer systems, second we perform simulations of the simple Doi-Edwards type reptation model. This 
model is one of the simplest models for entangled polymers, yet describes some essential features of entangled polymer 
dynamics qualitatively. The details of the model and the numerical scheme used in the simulations are described in 
Appendix [B] Although several important mechanisms (such as the constraint release) [!rt3j are not considered in 
this model, we believe that the incorporation of these mechanisms do not change the RF behavior of the TAMSD 
qualitatively. Therefore, if it gives qualitatively the same result as the slip-spring model, we consider that the RFA 
can be applied to most of models for entangled polymers. 

For simplicity, we employ the discrete version of the tube model. (In the followings we call this model as the 
discrete reptation model.) As described in Appendix IBl we model the tube as a freely jointed chain type model. A 
tube consists of Z tube segments and each tube segment has a fixed length a. (The number of segments in a tube is 
fixed to be constant and thus the contour length fluctuation (CLF)[l[ is absent in this model.) The polymer chain 
inside the tube moves (reptates) forth or back randomly, and we express the characteristic time scale of this random 
motion as tj = Qa? /UbT (Q is the friction coefficient of a segment in the logitudinal direction along a tube). We use a 
and T\ as elementary length and time scales. The time evolution of a tube is modeled as a random jump type process. 
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The chain inside the tube reptates forth or back by a given jump rate, and then a segment at an end is destructed 
and a new segment is constructed at another end. We perform kinetic Monte Carlo simulations for various different 
values of Z . (Z is proportional to the number of chain segments A.) The number of polymer chains (or tubes) is 
1024 for each simulation, and the initial state is sampled from an equilibrium distribution. Because this model lacks 
the CLF mechanism, the MSD of the CM becomes the following form[23j], instead of Eq. (|2|): 

([rcu(t) -r CM (0)] 2 } ext. (6) 

The EAMSDs of CMs are shown in Fig. All EAMSDs are proportional to t in all the time regions, which is 
consistent with Eq. ([6]). Therefore, in the discrete reptation model, we cannot obtain any characteristic time scales 
from the EAMSD at all. 

Interestingly, even in this case, the RE of the TAMSDs shows a clear crossover. As in the case of the slip-spring 
model, an appropriate value should be chosen for A in the REA. Here we set A/t; = 10. (We show the effect of A 
value in Appendix IAT) As shown in Fig. [5] the RF is fitted well with Eq. ([5]). Unlike the case of the slip-spring model, 
the exponent a is almost zero, which is different from the results of the slip-spring model. Moreover, we observe that 
the crossover time strongly depends on Z in a similar way to the case of the slip-spring model. The crossover time 
t c obtained from the RF of the TAMSDs are shown in Fig. [6] The obtained r c data almost coincide with the longest 
relaxation time Td or the reptation time r rep = t;Z 3 /7t 2 . This becomes clearer if we rescale the time t by the reptation 
time Trop- Figure [7] shows the RFs of TAMSDs (the same data as Figure [5|) for the rescaled time t/r re p. Interestingly, 
the rescaled RF data for different Z collapse into one master curve if Z is sufficiently large, except the short time 
region. 



III. DISCUSSIONS 



Both the slip-spring model and the discrete reptation model in Section |TT] give qualitatively the same crossover 
behavior. However, there are some differences between the REA results of the slip-spring and discrete reptation 
models. For example, the silp-spring model exhibits the power-law type behavior of R(t; A) for t < r c while the 
discrete reptation model exhibits almost constant R(t;A) for t < r c . We expect that the differences reflect the 
detailed relaxation mechanisms of models. (In the case of the CTRW, R(t; A) reflects some informations of traps 16].) 

The discrete reptation model is the simple model and the reptation is only the relaxation mechanism in the model. 
On the other hand, the slip-spring model has other relaxation mechanisms such as the Rouse motion of subchains 
and the contour length fluctuation. The incorporation of additional relaxation mechanisms other than the reptation 
modulates several dynamical behavior [H-[3l Izij . In the reptation model, the dynamic equation of the CM is described 
by using the end-to-end vector[23j]. Thus the RFs can be also modulated by various relaxation mechanisms, for 
example, through the dynamics of the end-to-end vector. If our expectation is correct, we will be able to extract the 
information about the relaxation mechanisms from the REA. However, currently we have no analytical theory for the 
RFs in entangled polymer sysmtes. Unfortunately, it seems to be difficult to obtain analytical expressions of RFAs of 
CMs, and the development of theoretical tools for RFs is desired. 

Although we have examined only two models and some mechanisms are not considered in them, we believe that we 
will observe similar behavior of TAMSDs for other models of entangled polymers. Judging from our simulation data, 
the reptation motion is the most important for the behavior of the RFs. Therefore, even if various mechanisms which 
affect the relaxation behavior are incorporated, we expect that the RFA still gives the information about the long 
time relaxation due to the reptation. It should be pointed that several mechanisms such as the CR affect the longest 
relaxation time Td [U, [H, [24[ • We expect that the crossover time r c will be affected in a similar way, and thus the RFA 
will give qualitatively the same results for more elaborated models. The numerical factor would be largely affected by 
mechanisms such as the CLF and the CR. The discrepancy between Td and r c observed in the slip-spring model may 
be due to such additional relaxation mechanisms. From these considerations, we expect that the RFA performed in 
Section |H] can be applied for other models, including the Kremer-Grest modelQ or highly coarse-grained models@,||j]. 

Here we shortly discuss the power-law exponent a of the slip-spring model (Eq. (0)). As we have shown in 
Section III Al it depends on A. As far as the authors know, there is no relaxation mechanism which gives such 
A-dependent power-law exponent. Therefore, it would be physically natural to interpret that this A-dependent 
power-law exponent is apparent. As we can observe in Fig. [2fb), the RF becomes proportional to i~ - 5 for a very 
small t region. The region where we have observed the power-law exponent a is the transient region. If N is relatively 
small, such a transient region is affected by the short time region, and the fitting may give apparent power-law 
exponent which is different from the true value. (This situation may be somewhat similar to the power-law for the 
longest relaxation time, Td <x N 3A , which is considered to be apparent behavior [lj].) 

One advantage of the RFA is that it can be performed only with the CM trajectory data. As we mentioned, for some 
cases, the CM trajectory data are more suitable for the calculation of the long time relaxation behavior. The accurate 
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computation of the stress autocorrelation functions in multi chain models is not easy. In highly coarse-grained models 
(such as the RaPiD model), the entangled polymers are modeled with very limited degrees of freedom, and some 
information such as the end-to-end vectors cannot be calculated. The RFA can be applied even for such cases. Of 
course, the RFA will not always be the optimal method, and other methods which utilize conventional quantities 
will be better for some cases. The RFA of the CM trajectories should be understood as a complementary method to 
conventional ones. Here we note that in conventional methods to determine the longest relaxation time th, dynamical 
quantities are utilized after the ensemble averages are taken. On the other hand, the RFA of the TAMSDs uses 
the time-averaged quantity before the ensemble average is taken. The RFA for other quantities, which are usually 
analyzed after the ensemble averages are taken, may be possible. 

Another advantage of the RFA will be that we can utilize it to validate some phenomenological coarse-grained 
models. As we will discuss in the next subsection, there are many possible ways to model the dynamics of entangled 
polymers. The RF behavior of TAMSDs will be useful to check the validity of a phenomenologically constructed 
dynamical model. That is, if the target model reasonably reproduces the dynamics of entangled polymers, it should 
exhibit the crossover behavior and the crossover time r c should be comparable to the longest relaxation time Td 
(estimated from other physical quantities). For such a validation, usually the ensemble-averaged quantities (such as 
the relaxation modulus) are utilized. We expect that the RFA, which utilizes a time-averaged quantity, may become 
a complementary validation tool. 

IV. CONCLUSIONS 

In this work we have shown evidences that the RFA of the TAMSDs of CMs of polymer chains can be utilized to 
study the long time relaxation behavior of entangled polymer systems. The information of the longest relaxation time 
(the disengagement time) Td is successfully extracted as the crossover time t c in the RF of the CM trajectories. 

The RFA can extract the characteristic long time relaxation behavior of entangled polymers for two different 
models (the slip-spring model and the discrete reptation model), although the behavior somewhat depends on the 
model details or parameters. For the discrete reptation model, the crossover time determined from the RF of the 
TAMSDs almost coincide with the longest relaxation time determined from the shear relaxation modulus. 

Although in this work we provided results only for two models, we believe that our method can be applied for 
other models of entangled polymers, including molecular dynamics models and highly coarse-grained models. Further 
analyses for various simulation data and development of theories for RFs in entangled polymer systems are demanding. 
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Appendix A: Effect of Time Difference A 

In this appendix, we examine the effect of the time difference size A on the RF of the TAMSDs, R(t; A). Figure [8] 
shows the RF for N = 80 and different A values in the slip-spring model. The forms of R(t; A) data depend on A 
rather strongly. If A is too large (such as A/tq = 1000), we cannot determine the crossover time r c from the RF 
data. This may be because the A value becomes close to and the crossover behavior may be smeared out. If A 
is too small (such as A/to = 1), the R(t; A) value becomes very small and thus it becomes difficult to determine r c 
accurately. Thus A is required to be an intermediate value. In this case, the curve for A/tq = 10 seems to be better 
than other curves. For other TV cases, we observe a similar trend. Therefore, in this work we employ A/to = 10 to 
determine the crossover time r c . Physically, this value is comparable to the entanglement time T e . (We note that even 
if we employ other value for A, for example A/to = 100, the result in the main text is qualitatively not affected.) 

Figure [3] shows the RF for Z = 40 and several different A values in the discrete reptation model. As clearly observed 
in Fig. IH1 the R(t; A) data for relatively small A (A/77 S 10) are almost the same. Thus we consider that there is a 
threshold for A/77 an d if A is smaller than the threshold, R(t; A) is insensitive to A. The threshold value depends on 
several parameters such as Z. To determine the crossover time t c for different Z values, we should choose A smaller 
than the threshold for all the Z values examined. In this work, we employ A/77 = 10, which is sufficiently small for 
all the examined parameters in the main text. In this case, 77 is comparable to r e and thus this A is slightly larger 
than r e . 
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In conclusion, the results in the main text are not so sensitive to the values of A, as long as A is selected to be an 
intermediate value. Results in this appendix implies that the A value should be selected comparable to (or slightly 
larger than) the entanglement time r e . Although the reason why such A values work well is not clear, this may be 
useful to estimate the optimal value of A when one performs the RFA for other simulation results. A qualitatively 
similar result is obtained by a theoretical analysis for the CTRW[16]. 



Appendix B: Discrete Reptation Model 



In this appendix, we provide the details of the discrete reptation model and the numerical scheme used for simu- 
lations. In the Doi-Edwards type reptation model[l[, the polymer chain is assumed to be confined in a tube which 
consists of Z discrete steps. Each step has a constant size a, and thus the tube can be regarded as a freely jointed 
chain like object[l[. In this work we express the positions of tube ends and kinks by {Ri} (i = 0,1,2, ... , Z). The 
chain inside the tube can move only to the direction along the tube (which we call the longitudinal direction) , and this 
reptation motion is modeled by the Langevin equation. However, the Langevin equation is not suitable to describe 
the dynamics of discrete objects, especially when we perform numerical simulations. In this work, we employ the 
discrete jump dynamics to model the reptation motion. 

In the discretized reptation model, the polymer chain can move forth or back along the chain by the constant 
segment size a with one step. We ignore the contour length fluctuation (CLF) effect and the constraint release (CR) 
effect. Thus the dynamics of the chain is described only by the reptation motion described here. The jump rates are 
proportional to Z^ 1 , because the jump rate corresponds to the longitudinal motion of a whole chain along the tube 
and the effective friction coefficient of a chain is proportional to Z^ 1 . The forth or back jump rates (transition rates) 
are given as 

W± = -±- (Bl) 

where 77 is the characteristic longitudinal diffusion time defined as 17 = (ia 2 /ksT is the friction coefficient of a 
segment in the longitudinal direction). Subscripts "+" and "— " represents the forth and back directions, respectively. 
This transition rate model successfully reproduces the diffusion coefficient of a chain along the tube (in the longitudinal 
direction), D t = k B T/ZQ^. 

To avoid ambiguities or problems associated with the dynamical mapping of Monte Carlo simulations (25i [26| . 
we employ the kinetic Monte Carlo scheme [26l - [28| which can handle the time evolution of discrete stochastic jump 
processes exactly. We start from the equilibrium tube conformation {Ri} at t — to. To generate random numbers, 
we employ the Mersenne twister pseudo random number generator [29j. The kinetic Monte Carlo scheme evolves the 
system from time t = tj to t = t 3+ i where j is the number of Monte Carlo steps. We need the total transition rate 
Wtot and the (normalized) probabilities of transitions P± , to simulate the time evolution. They are given as 

W to t =W+ + W- = (B2) 

The direction (forth or back) is randomly selected according to the probability P±. From Eq. (IB3|) the direction is 
selected with the same probability, and this can be easily realized by using a uniformly distributed random number. 
The time step size is sampled from the exponential distribution. The time at the (j + l)-th step becomes 

*j+i = *i-ri7-hu (B4) 

vvtot 

where u is the random number sampled from the uniform distribution in (0, 1). The positions of tube ends and kinks 
at the (j + I)-th step are given as 

n( . fR t+ i(t 3 ) (i<Z) JfloteJ + on (i = 0) 

R ^^-\ Rz[t]) + an {l = Z) (i>0) (B5) 

where n is the random vector on the three dimensional unit sphere (\n\ = I). The time series of the tube conformation 
are obtained by successively iterating the kinetic Monte Carlo time steps. 
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The segments of the polymer chain are assumed to be uniformly distributed along the tube. Then the CM position 
is calculated from {Ri} as 

z=0 

For the discrete reptation model, we define the longest relaxation time Td by using the shear relaxation modulus 
G(t) at the long time limit, 

-t'-H-^ <«-«»■ < B7 > 

(This method to determine Td is consistent with one used for the slip-spring model.) We calculate the shear relaxation 
modulus by the linear response formula: 

G(t) =c o k B T(a xy {t)a xy {0)). (B8) 

Here cq is the spatial average chain number density and a xy is the xy component of the stress tensor. The stress 
tensor cr is calculated as 

^ = ^£ yVJli+i - Ri)(Ri+i - Ri) - (Z + l)k B Tl. (B9) 

a u 

We note that Eqs. (|B7[) - (|B9p give the shear relaxation modulus which is proportional to the tube surviving 
probability [lj]. Therefore it is almost obvious for our model that the longest relaxation time Td coincides to the 
reptation time r rcp = QZ 3 a 2 /-K 2 k B T — Z 3 ti/tt 2 . 
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Figure Captions 

Figure Q] Ensemble-averaged EAMSDs of CMs in the slip-spring model for N = 10, 20, 40, 80, and 160. Dashed lines 
represent curves which are proportional to A 1 and A 1 / 2 , b and To are the size and characteristic time of a segment, 
respectively. 

Figured Relative fluctuations of TAMSDs of CMs in the slip-spring model for A/r = 10 and N = 10, 20, 40, 80 and 
160, where To is the characteristic time of a segment, (a) N = 10 and 20, and (b) N = 40, 80, and 160. Dashed line 
represents a curve proportional to t^ 1 / 2 . The power law exponents for the short time regions are a — 0.31, 0.24, and 
0.19 for N = 40, 80, and 160, respectively. 

Figure [3j The longest relaxation time Td and the crossover time t c in the slip-spring model. Td is determined from the 
shear relaxation modulus data [22j whereas t c is determined from the RF data of TAMSDs for A/tq — 10. Dashed 
lines represent fitting curves for large N (t<2 oc N 3AS and t c oc TV 3 51 ). 

Figure SI Ensemble-averaged MSDs of CMs in the discrete reptation model for Z = 10,20,40,80, and 160. a is the 
step size of a tube segment and t/ is the characteristic time of the longitudinal motion of a segment along the tube. 

Figure [5] Relative fluctuations of TAMSDs of CMs in the discrete reptation model for A/t; = 10 (see Appendix |A"|) 
and Z — 10, 20, 40, 80 and 160. The dashed line represents a curve proportional to t~ x / 2 . t\ is the characteristic time 
of the longitudinal motion of a segment along the tube. 

Figure IH1 The longest relaxation time t^ and crossover time t c in the discrete reptation model. Td is determined from 
the shear relaxation modulus and t c is determined from the RF data of TAMSDs in the same way of Fig. [3] The 
dashed line represents the reptation time t iop /ti — Z 3 /it 2 . 

Figure [7J Rescaled RFs of TAMSDs of CMs in the discrete reptation model for different values of Z. The data are 
the same as Figure [5] but the observation time t is rescaled by the reptation time T rcp = Z 3 ti/tt 2 . All the data points 
collapse into one master curve except the short time region or the small Z data (Z — 10, in this case). 

Figure [S] Relative fluctuations of TAMSDs of CMs in the slip-spring model for N = 80. The dashed line represents 
a curve proportional to i" 1 / 2 . 

Figure [S] Relative fluctuations of TAMSDs of CMs in the discrete reptation model for Z — 40. The dashed line 
represents a curve proportional to t~ x / 2 . 
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